Wigner Crystallization in Electron-Hole Double Layers: 
Evolution from Long- to Short-Range Forces. 
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We study fermions in two dimensions interacting via a long range 1/r potential for small particle 
separations and a short range 1/r 3 potential for larger separations on length scales £. We computed 
the energy of the Wigner crystal and a homogeneous Fermi liquid phases in a variational approach 
and determined the phase diagram as a function of density and £ at zero temperature. We discuss 
the collective modes in the FL phase, finite temperature effects on the phase diagram, and possible 
experimental realizations of this model. 
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I. INTRODUCTION AND RESULTS 



Inhomogeneity plays an important role in many highly 
correlated materials 1-8 . For example stripe phases 
which break translational symmetry of real space in 
one direction have been observed in high temperature 
superconductors 4 . Similarly a Wigner crystal (WC) 
phase 1 which has only discrete translational symmetry 
of triangular lattices and six-fold rotational symmetry 
has been observed in electrons on the surface of Helium 2 
and in ultra clean two dimensional hole gases 3 . 

For interactions which decay as a power law l/r Q the 
potential energy scales with density as n a / 2 in two di- 
mensions (2D) and the kinetic energy as n. The ratio 
of the potential to the kinetic energy scales as rt Q / 2-1 
which becomes rt -1 / 2 for Coulomb interactions, a con- 
stant for 1/r 2 and n 1 / 2 for 1/r 3 potentials. A stable WC 
crystal phase occurs when interactions are of the order 
of kinetic energy 1 . Hence a WC is expected in the high 
density regime or low density regime depending on the 
power law decay of interactions. For Coulomb forces this 
suggests that there is a phase transition at decreasing 
densities from a homogeneous Fermi liquid (FL) phase, 
which is conducting, to a WC phase, which is insulating. 
In this regime, perturbative methods fail thus one must 
resort to numerical studies 5,6 . 

In reality, many important systems do not have pure 
power law interactions. In two dimensional electron gases 
(2DEGs) in inversion layers or semiconductor quantum 
wells, the presence of a nearby gate modifies the Coulomb 
interaction. Image charges will form in the gates and 
they screen the interactions such that at short distances 
it is 1 jr and at long distance is 1/r 3 with a length scale 
given by the distance to the gate. The effects of such 
long-range forces dramatically changes the energetics of 
the 2DEG. A phase diagram was estimated in Ref. 7,8 . 
Recently the effects of long-range forces on superconduc- 
tivity have been explored too 9 . A similar experimental 
scenario occurs in ionic liquid transistors 10 where an elec- 
trolyte is used as a dielectric in a standard field effect 
transistor configuration. Positive and negative charges 
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FIG. 1. (Color online) Phase diagram of spinless fermions in 
two dimensions with 43 particles showing the phase bound- 
ary between the Wigner crystal and uniform FL phases as 
function of the screening length £ in units of ao and density 
n = \j (jia\r 2 s ). Also shown is the £ = r s line which divides 
the Coulomb 1/r regime (right) from the 1/r 3 regime (left). 
The dashed (red) lines represent the analytical asymptotes of 
the boundaries (see text). 



accumulate at opposite ends of the electrolyte. Finally, 
recent progress in cooling techniques have allowed the 
study of degenerate dipole gases 11-14 . Such objects in- 
teract as 1/r 3 at large separations. At short distances 
the interaction is modified by the interactions between 
the electronic clouds. The common feature of these sys- 
tems is that there is a length scale in the interactions. 
They are composed of charges of opposite sign which are 
located in layers separated by a dielectric. As such, the 
interparticle potential crosses over between two different 
power law regimes. 

In this work, we study a minimal model in a varia- 
tional approach of spinless fermions interacting with a 
1/r (1/r 3 ) potential at short (long) distances in two di- 
mensions, see Eq. 2. Our results agree with Ref. 7,8 in 
the study of 2DEG with a nearby gate, but we provide a 
microscopic model Hamiltonian from which the phase di- 
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agram arises. Moreover we find that this phase diagram 
should qualitatively describe other systems with similar 
crossover. 

The corresponding zero-temperature (T = 0) phase 
diagram is presented in Fig. 1 as a function of the sepa- 
ration £ and density n parametrized by r s = 1/ (do y/im) ■ 
The FL and WC phases are indicated in the phase dia- 
gram and labeled according to the dominant interaction 
regimes. For fixed £ > £ c there is a regime of densities 
where the Wigner crystal is energetically more favorable 
than a uniform FL. We also show in dashed(red) lines 
the asymptotic forms of the phase boundaries from our 
numerical simulations. In the regime dominated by l/r 3 

the asymptotic phase boundary corresponds to £ ~ r~J 2 \ 
while in the Coulomb case it is r s — constant (~ 13). In- 
tuitively, for ~ 1 /r 3 potentials (left to the £ = r s line) the 
kinetic (~ l/r 2 ) and potential energy ~ e 2 £ 2 /rf are com- 

1 /2 

parable when £ ~ r s , whereas for the l/r regime, the 
potential energy scales as ~ e 2 /r s and these energies are 
comparable when l/r 2 is a constant (pa 5.9 x 10 -3 ). No- 
tice that there is a maximum value of 1 /£ 2 , corresponding 
to the condition £ c pa 15oo, below which the WC is more 
favorable for a wide range of r s . The existence of the 
maximum is a general feature that follows also from the 
above physical argument; if £ is too small the potential 
energy cannot be simultaneously of the same order as the 
kinetic energy. The limit £ — > oo is not directly acces- 
sible with our method. Extrapolating we find that the 
FL phase is more stable than the WC phase for r s < 13. 
For l/£ 2 below the critical maximum and increasing den- 
sity from zero to finite values, we find a first order phase 
transition from a 'dipolar' FL to a WC. With additional 
increases of the density beyond the £ = r s line, the WC 
crosses over to a Coulomb WC. Upon further increase 
of the density, the WC phase transitions to a FL. This 
reentrant behavior for the FL phase is a remarkable char- 
acteristic of any 2D system interacting with an effective 
potential which crosses over from long to short-range. 
Fig. 1 constitutes one of our main results. 

In what follows we show how this phase diagram is 
obtained. In section II we explicitly define the Hamilto- 
nian and in section III we explicitly compute the energy 
of the homogeneous FL and discuss its collective modes 
as they are modified by the presence of £. In section IV 
we compute the energy in WC phase with a variational 
wave function. In section V we compare the energies of 
these two states in different density regimes and deter- 
mine the ground state of the system. A discussion of the 
scaling of the WC energy with density is given along with 
a schematic finite temperature phase diagram. 

II. MODEL 

We start from a microscopic Hamiltonian 
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FIG. 2. (Color online) Energy per particle vs density (in units 
of particles per £ 2 ) of the WC and FL phases. Here £ = 30ao. 
Inset: E wc (s) vs size of the wave function. 

where the first term corresponds to the kinetic energy 
and the second term 

. . 2e 2 2e 2 

V r = / 2— 72W2 . 2 

describes the interactions of, say, electrons and their im- 
age charges located at the gate. The variable £ is the 
'screening' length and we treat it here as a phenomeno- 
logical parameter. Physically it also could represent the 
classical effective dipole size. In the first case one would 
include the energy of the capacitor energy that forms be- 
tween the gate and the 2DEG. However such effects are 
not considered here as they likely give small corrections 
to the mean field diagram 7 ' 8 , e is the dielectric constant 
of the medium. Note that for fixed density in the limit 
£ — > the interaction becomes V(r) = 2e 2 /er twice the 
Coulomb potential whereas in the limit £ — > it becomes 
V(r) = 2(e£) 2 /er 3 to lowest non-vanishing order in £. 

For both FL and WC phases, we consider a many-body 
ground state wave function to be 

#(n,r 2) ...,riv) =-4.[<Mri)---<Mrjv)], ( 3 ) 

where ./V is the number of particles, and A in a oper- 
ator that describes the anti symmetrization necessary 
for indistinguishable fermions. Results for indistinguish- 
able bosons and for classical distinguishable particles are 
straight forward extensions. For simplicity, we assume 
that the spin degrees of freedom are frozen, such that the 
wave function describing individual particles 4>i{yi) are 
either plane waves in the FL phase, or Gaussian centered 
at individual sites of a triangular lattice in the Wigner 
crystal phase. 

III. FERMI LIQUID PHASE AND COLLECTIVE 
MODES 

To begin, let us consider the homogeneous FL phase. 
Our system has two length scales. The first scale is £ 
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while the second is the interparticlc distance kp 1 set by 
the density n = kp/ (47r) , where kp is the Fermi momen- 
tum of the FL with Fermi energy Ep = kp /2m. The rela- 
tive size of these length scales determines the range of the 
potential. For £/ (anr s ) <C 1 particles interact with a 1/r 3 
potential, whereas for £/(aor s ) 3> 1 via the Coulomb 1/r 
potential, see Fig. 1. Explicitly, the energy per particle 
of a uniform polarized uniform FL is 



interaction of the form Eq. 2 we obtain, 

,2\ 1/2 



E FL _ 2 | 4£ 2 
r 2 £ 

, 8£ 



Io(4£/r,)-I a (4£/r a ) 



l-io^/r.J+La^/r.) 



(4) 



where £ = C/ a o is scaled by the Bohr radius ao = 
eh 2 /me 2 , and the energy scale E = e 2 /2ao- The func- 
tions I n [x) are the n th order modified Bessel function 
of the 1st kind and L n (x) are the n th order Struve func- 
tions. In obtaining Eq. (4), we used the fact that the two- 
body interactions shown in Eq. (2) have Fourier trans- 
form V(q) = (47re 2 /e<7)(l - e"^). 

The first term in Eq. (4) is the kinetic energy and 
scales as the density n. The second term in Eq. (4) is the 
Hartree contribution, which also scales as n. In the limit 
£ -> co, E FL /E Q = 2/r 2 + 4£/r 2 - 32/(37rr s ) is formally 
divergent as interactions are long-range 1/r. In a 2DEG 
the positive background charge cancels the Hartree term 
and in this case, Epl/Eq = 2/r 2 — 32/(37rr s ), in agree- 
ment with known results 15 . The interaction term is twice 
that of a 2DEG, see our Eq. (2). The last term within 
the square brackets is the Fock contribution which be- 
comes — 4^/r 2 + 512£ 2 /(457rr 3 ) as £ — > and leads to 
Epl/Eo = 2/r 2 + 512£ 2 /(457rr 3 ) + 0((?) in this 1/r 3 
regime. To 0(£) the Fock term cancels the Hartree term 
for pure 1/r 3 potentials. Physically, the system behaves 
as charge neutral FL of dipoles. We note that including 
a background to enforce charge neutrality as in 2DEGs 
with 1 /r potentials leads to negative energies and a self- 
bound system. In our case however, the energies are pos- 
itive and the system has positive pressure. The energy 
of the uniform FL is shown in Fig. 2 for £ = 30ao. 



A. Collective modes 

The length scale £ in the potential introduces impor- 
tant modifications to the collective excitations in the FL 
phase. The Fourier transform of Eq. (2) at zero mo- 
mentum is well-defined, V(0) = 47re 2 £. This means that 
at sufficiently long distances the interaction is effectively 
short range 1/r 3 . This introduces important modifica- 
tions to the energetics of the system 16 . For example, the 
Hartree and Fock terms must be included in the self en- 
ergy to have a conserving approximation to the Green 
functions 17 . The collective modes are given by the zeros 
of the dielectric function e(q,u>) = 1 — V(q)H(q,uj). A 
non-conserving RPA calculation 18 gives a good approx- 
imation to the behavior of the collective modes. For an 
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This result is valid at weak interactions and small wave 
vectors. As we see, the system supports zero sound 
modes whose group velocity is proportional ~ £}/ 2 . The 
speed of the zero sound diverges as £ — > oo which means 
that the power expansion in q/fl breaks down. Indeed, 
the plasmon dispersion relation has been extensively 
studied in 2DEGs in GaAs 19,20 semiconductor quantum 
wells. The Fourier transform of the pure Coulomb po- 
tential V(q) ~ 1/q diverges at zero momentum and leads 
to gapless ~ Jq mode dispersion 20 . 

According to the variational principle, the energy ob- 
tained using the many-body wave function described in 
Eq. (3) is a rigorous upper bound for the true ground 
state energy of the system. In the following analysis, we 
perform a variational calculation of the ground state en- 
ergy of the Wigner crystal phase and compare it with 
that of the Fermi liquid phase obtained in Eq. (4). 



IV. WIGNER CRYSTAL PHASE 

In the Wigner crystal phase, the single particle wave 
functions are localized at sites i of a 2D triangular lat- 
tice, 0j(r) = (l/s^/Trjcxp [— (r - R ?; ) 2 /(2s 2 )] , where Rj 
is the site position and s parametrizes the 2D size of 
the wave function. These single particle wave functions 
are approximately orthonormal, since the overlap at dif- 
ferent sites is exponentially small, J d 2 r cf>*(r)(j)j(r) = 
cxp [-(Hi - Hj) 2 /4s 2 ] = exp [-i??-/(4s 2 )] , provided 
that the separation R^ between sites i and j is much 
larger than the extent s of the single particle wave func- 
tion, that is Rij 3> s. We denote the lattice spacing 
by I and consider the regime of weakly ovelaping single 
particle wave functions, where s < 1/2. The localized na- 
ture of the many-body wave function is reasonable only 
in this regime. Explicitly, the energy per particle in the 
WC phase is, 

En S l In 



(W^ 2 Erfc(^)J F 2 ( S "), (6) 

where the functions Fi(y) are lattice sums given by 
-^i(^) = 2~2i^j Jo(kRij)/N appearing in the second term 

of Eq. (6), and F 2 (s) = £ ¥j ,e~ R ^' 2 ^ /N appearing 
in the third-term of Eq. (6). s — s/clq, k — kan, 
Rij = Rij /ao and J„ is the n th order Bessel function 
of the first kind. The first term is the zero point motion 
due to the localization of the particles at the lattice sites. 
This term favors large single particle wave functions. The 
second and third terms are the Hartree-Fock (HF) con- 
tribution which is non-monotonic in s. The extension to 
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particles with Bosonic statistics results in a change of the 
sign of the Fock term whereas for classical particles only 
the Hartree term is present. 

We comment on the limiting behavior of Eq. 6. For 
£ — > oo the system interacts with a Coulomb potential. 
In this limit, 



/^Coulomb 
wc 

En 
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£ e -^ 2 /o(^/4s 2 ) 



sN 



(7) 



In the regime of non-overlapping wave functions l/2s^$> 
1, the asymptotic form of the Bcsscl function In(x) ~ 
e x j\j2nx. As we see, the series in the second (Hartree) 
term diverges when N — > oo as expected since we have 
not included a neutralizing background charge. In this 
limit, taking a finite N gives a well defined energy which 
is non-monotonic in s. We find that this behavior extends 
to the £ 7^ oo regime and produces a WC state energy 
Evrd 3 ) that does not have a minimum as a function of 
s when the density is either too small or too large. On 
the other hand, in the limit of £ — > we obtain that the 
0(£) contribution of the Fock term cancels the Hartree 
term (as in the FL phase) and the first non-vanishing 
term is 0(£ 2 ). Since the resulting expression is not very 
illuminating, we omit it here. 



V. DISCUSSION AND CONCLUSIONS 

For £ ^ oo, the integrals appearing in the analytical 
expression of Eq. (6) need to be calculated numerically. 
We have computed numerically the energy per particle as 
a function of the density n = 2/(Pv / 3) for the WC phase 
in a 2D triangular array of 43 and 55 particles. In Fig. (1) 
we show the resulting phase diagram for 43 particles, but 
wc would like to emphasize that the qualitative behavior 
found is essentially unchanged with 55 particles. How- 
ever, the calculation of the phase boundaries for the full 
range of parameters £ and r s is very intensive in the limit 
of large number of particles. 

We establish the phase boundary between the solid 
and liquid phases by minimizing the energy in the solid 
phase with respect to the variational parameter s. A 
typical variation of the energy as a function s is shown 
in Fig. (2). When there is a local minimum of E wc (s) 
located at s = sq, we compare the energy E wc (sn) with 
the energy of the uniform screened FL with the same 
density. We find that the energy of the WC is always 
lower than the FL energy in the regime of tested values 
of £. However a local minimum of E wc (s) does not exist 
for densities below a minimum n c ,min(£) and for densi- 
ties above a maximum n c ,max{£), where the WC phase is 
unstable. This establishes that the WC phase is unsta- 
ble for densities n satisfying the condition n < n Ctm i n {£/) 
and n > n c , mal ((), in which case, the FL phase is the 




FIG. 3. (Color online) Scaling of the energy of the Wigner 
crystal (WC) near n Ct min and n c , max for 43 particles. Here 
£ = 30a . 



stable one. In the limit £ — > oo and fixed density, we 
find that the critical value of density for the solid-liquid 
phase transition occurs when r s = 13 in the case of 43 
and changes little for 55 particles. 

A phase transition from a ferromagnetic FL to a ferro- 
magnetic WC phase in a 2DEG is expected to occur for 
r s = 29 according to recent Monte-Carlo simulations 6 . 
In the £ — > oo limit we find such a phase transition for 
r s = 13, but our interaction is twice that of pure Coulomb 
with finite size samples. We emphasize that our results 
does not apply in this regime as interactions are of pure 
Coulomb character and a neutralizing background charge 
must be explicitly considered along with numerical meth- 
ods. Here we are interested in the regime of £ < oo which 
we expect to be qualitatively correct. 

It is interesting to note that in our simulations with 
43 and confirmed with 55 particles the optimized energy 
values in the Wigner crystal phase can be fitted with a 
power law 



E, 



,4/3 



(8) 



for values of density close to the boundary with the FL 
phase, n — > nf :min in the regime of predominantly 1/r 3 
forces. See also Fig. 1. This is shown explicitly in 
Fig. 3 for £ = 30ao- Additionally, in the tested range 
of 20 < £/ao < 40 we find that the exponent is 1.33 
to within 1% error. Similar scaling behavior has been 
studied in 2DEGs with a screened 1/r 3 potential and 
spin-orbit interaction 21 . Similarly the optimized energy 
in the Wigner crystal can be fitted by a linear power 
E wc = C\n — n C: , nax \ + D for densities close to the FL 



phase n 



,., where C and D are constants but in 



the 1/r regime. In the tested range of 20 < £/ao < 40, 
wc find that the exponent is 1.00 to within 1% error. 
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FIG. 4. Schematic variational phase diagram of the electron- 
hole double layer as function of density ~ 1/rf, temperature 
T and layer separation £ in appropriate units. V is the ratio 
of the potential to kinetic energy. See text for details. 



Finite temperature mean field phase diagram of 
a screened 2DEG 



Here we estimate the mean field phase diagram cor- 
responding to a 2DEG in the presence of a nearby gate 
at finite temperature. As discussed above we assume 
this system can be modeled as an electron-hole double 
layer with separation £. The phase diagram in the l/£ 2 - 
density plane had been estimated before 7 from general 
considerations. At finite temperatures, we expect the 
melting of the WC to a liquid-like state (via Hcxatic or 
bubble phases 8 ). The limit £ — > oo is not directly acces- 
sible with our method but can be estimated from general 
considerations 20 . In Fig. 4 we show the phase diagram in 
the (T, n, l/£ 2 ) parameter space. In the l/£ 2 = plane 
the thermodynamics of the classical 2DEG is determined 
by the quantity T which is the ratio of the potential to 
the kinetic energy per particle. For a classical 2DEG this 



ratio is T = (7rn) 1 / 2 e 2 /T. For T < 1 the kinetic energy 
dominates and the system behaves as a classical liquid. 
For r> 1 the Coulomb potential dominates and we ex- 
pect a classical WC. The phase boundary is defined by 
T = T c - 100, see Ref. 2 . 

For temperatures below Ep ~ n ~ 1/r 2 the Pauli 
principle and quantum correlations drive T down with 
increasing density 20 , e.g., along the 1/r 2 line in Fig. 4, 
r ~ ji -1 / 2 in the Hartree-Fock approximation. This de- 
fines a critical density above which T < 1 and a quantum 
fluid is recovered. The full melting line and the quan- 
tum/classical crossover in the plane l/£ 2 = is sketched 
in Fig. 4. Connecting the WC phases of the l/£ 2 = and 
T = planes is a dome where a WC-like phase should be 
stable. Within this dome there exists at least 4 distin- 
guishable phases corresponding to the quantum vs clas- 
sical and 1/r vs 1/r 3 regimes. Outside the dome there 
are only quantum and classical liquid- like phases. 

In conclusion, we have obtained from a microscopic 
model, the phase diagram of fermions with an interac- 
tion which is smoothly deformed between two power law 
regimes which lie in the long-range and short-range cate- 
gory in 2D. We argue that such phase diagram is generic 
to any system with such crossover, e.g, 2DEG in semicon- 
ductor quantum wells with gates 7 , and polarized dipole 
gases when the size of the dipoles is taken into account. 
Our results have relevance to experiments on polarized 
dipole gases in 2D where a WC is favorable at high den- 
sities. We suggest that at such high densities one likely 
should consider corrections to the 1 /r 3 form of the po- 
tential. 



ACKNOWLEDGMENTS 

We thank the JQI-PFC (BMF) and ARO (W911NF- 
09-1-0220) (CSdM) for support. 



1 E. Wigner, Phys. Rev. 46, 1002 (1934). 

2 C. C. Grimes and G. Adams, Phys. Rev. Lett. 42, 795798 
(1979); E. Rousseau, D. Ponarin, L. Hristako, O. Avenel, 
E. Varoquaux, and Y. Mukharsky, Phys. Rev. B 79, 
045406 (2009). 

3 J. Yoon, C. C. Li, D. Shahar, D. C. Tsui, and M. Shayegan, 
Phys. Rev. Lett. 82, 1744 (1999). 

4 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75, 1201 (2003). 

5 B. Tanatar and D. M. Ceperly, Phys. Rev. B 39, 5005 
(1989). 

6 N. D. Drummond and R. J. Needs, Phys. Rev. Lett. 102, 
126402 (2009). 

7 B. Spivak, Phys. Rev. B 67, 125205 (2003). 

8 B. Spivak and S. A. Kivelson, Phys. Rev. B 70, 155114 
(2004). 



9 S. Raghu, E. Berg, A. V. Chubukov, and S. A. Kivelson, 
ArXiv:1111.2982vl. 

10 R. E. G. Ill and M. D. Sherrill, Phys. Rev. Lett. 5, 248250 
(1960). 

11 M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 
108, 215301 (2012). 

12 M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, 
Chem. Rev. 112 (9), 50125061 (2012). 

13 K. Sun, C. Wu, and S. Das Sarma, Physical Review B 82, 
075105 (2010). 

14 C. Lin, E. Zhao, and W. V. Liu, Phys. Rev. B 81, 045115 
(2010). 

15 A. K. Rajagopal and J. Kimball, Phys. Rev. B 15, 2019 
(1977). 

16 J. P. Kestner and S. D. Sarma, Phys. Rev. A 82, 033608 
(2010). 

17 B. M. Fregoso and E. Fradkin, Phys. Rev. B 81, 214443 



6 



(2010). 

L. P. Kadanoff and G. Baym, Quantum Statistical Mechan- 
ics (Benjamin, New York, 1962). 

M. Eriksson, A. Pinczuk, B. Dennis, C. Hirjibehedin, S. Si- 



mon, L. Pfeiffer, and K. West, Physica E 6, 165 (2000). 

20 E. H. Hwang and S. D. Sarma, Phys. Rev. B 75, 205418 
(2007). 

21 E. Berg, M. S. Rudner, and S. A. Kivelson, Phys. Rev. B 
85, 035116 (2012). 



